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to the cosmological lithium 
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The measurement of the cosmic microwave back- 
ground has strongly constrained the cosmologi- 
cal parameters of the Universe 1 . When the mea- 
sured density of baryons (ordinary matter) is com- 
bined with standard Big Bang nucleosynthesis 
calculations 2 ' 3 , the amounts of hydrogen, helium 
and lithium produced shortly after the Big Bang 
can be predicted with unprecedented precision 1 ' 4 . 
The predicted primordial lithium abundance is a 
factor of two to three higher than the value mea- 
sured in the atmospheres of old stars 5 ' 6 . With 
estimated errors of 10 to 25 % , this cosmological 
lithium discrepancy seriously challenges our un- 
derstanding of stellar physics, Big Bang nucleosyn- 
thesis or both. Certain modifications to nucleosyn- 
thesis have been proposed 7 , but found experimen- 
tally not to be viable 8 . Diffusion theory, how- 
ever, predicts atmospheric abundances of stars to 
vary with time 9 , which offers a possible explana- 
tion of the discrepancy. Here we report spectro- 
scopic observations of stars in the metal-poor glob- 
ular cluster NGC 6397 that reveal trends of atmo- 
spheric abundance with evolutionary stage for var- 
ious elements. These element-specific trends are 
reproduced by stellar-evolution models with diffu- 
sion and turbulent mixing 10 . We thus conclude 
that diffusion is predominantly responsible for the 
low apparent stellar lithium abundance in the at- 
mospheres of old stars by transporting the lithium 
deep into the star. 

Diffusive processes altering the elemental compo- 
sition in stars have been studied for decades 9 ' 11 . 
Evidence for their importance comes from 
helioseismology 12 and the study of hot stars with 
peculiar abundance patterns 13 . Among solar-type 
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stars, the effects of diffusion are expected to be more 
pronounced in old, very metal-poor stars. Given their 
greater age, diffusion has had more time to produce 
sizeable effects than in younger stars like the Sun. 
Detailed element-by-element predictions from models 
including effects of atomic diffusion and radiative 
accelerations became available a few years ago 14 , 
but these early models produced strong abundance 
trends that are incompatible with measurements of, 
in particular, the abundance of lithium common 
among stars of the Galactic halo over a wide range of 
metallicities (the so-called Spite plateau of lithium). 
However, the recent inclusion of turbulent mixing 10 
brings model predictions into better agreement with 
observations. 

According to the predictions from such model cal- 
culations, stars leaving the main-sequence (turn-off 
stars) are expected to show the largest variations 
relative to the composition of the gas from which 
the stars originated. Giant stars, however, have deep 
surface convection zones which erase most effects of 
diffusion and restore the original composition. One 
notable exception is lithium which disintegrates in 
layers with T > 2.1 million K. The destruction of 
lithium inside the star leads to a successive dilution of 
the surface lithium because the convective envelope 
expands when the star becomes a red giant. We 
performed spectroscopic observations specifically to 
test these model predictions. 

Globular-cluster stars have identical age and initial 
heavy-element composition, and thus measured 
atmospheric abundance trends with evolutionary stage 
are a signature of diffusion. We observed 18 stars 
along the evolutionary sequence of NGC 6397 with 
the multi-object spectrograph FLAMES-UVES on 
the 8.2-m Very Large Telescope in Chile - five stars 
close to the main-sequence turn-off point (TOP), two 
stars in the middle of the subgiant branch (SGB), 
five stars at the base of the red-giant branch (bRGB) 
and six red giants (RGB) - that represent specific 
stages in stellar evolution (Fig. 2). Even on 8-m-class 
telescopes, these observations are still challenging, 
requiring total integration times of 2-12 hours per star 
to obtain homogeneous data with high signal-to-noise 
ratios at high spectral resolution. 
These observational challenges explain in part why 
there have been surprisingly few studies looking for 
abundance trends between unevolved and evolved 
stars in globular clusters. One study 15 found lower 
abundances for subgiants in the very metal-poor 
globular cluster M92 than accepted literature values 



2 



Nature 442, 657-659 (10 August 2006) 



Mean stellar parameters, iron abundances and photometric quantities of the four groups of stars 



Group 


No. of star 


Tcff (K) 


log \q (cm s 


~ 2 )] 

/J 


logje (Fe)l 

oL V ' J 


£ (kms^ 1 ) 


TOP 


5 


6254 


3.89 




5.23 ± 0.04 


2.00 


SGB 


2 


5805 


3.58 




5.27 ± 0.05 


1.75 


bRGB 


5 


5456 


3.37 




5.33 ±0.03 


1.73 


RGB 


6 


5130 


2.56 




5.39 ± 0.02 


1.60 


Sun 


1 


5777 


4.44 




7.51 


1.00 




A T off (TOP -RGB) (K) 


A log g (TOP- 


RGB) 


Aloge(Fe) (TOP -RGB) 


Spectroscopy 


1124 




1.33 




0.16 ±0.05 


Stromgren, (v — y) 


1108 




1.38 








Broad-band, (V - I) 


1070 




1.38 









Table 1: Mean stellar parameters of the four groups of stars in the following evolutionary stages: TOP, SGB, bRGB and 
RGB (see text). The FLAMES-UVES spectra cover the spectral range from 4800 - 6800 A, have R = A/AA = 48 000 (where 
A is wavelength) and signal-to-noise ratios in excess of 80: 1 per pixel. The analyses are fully spectroscopic and line-by-line 
differential to the Sun. Typical errors on T e g, logp and the microturbulence £ for individual stars are, respectively, 150 K, 
0.15 and 0.2 km s _1 . The errors in log[e(Fe)] = log( Nf c /Nh) +12 are the combined values of the line-to-line scatter of Fe I 
and Fe II for the individual stars propagated into the mean value for the group. Between 20 and 40 Fe I and Fe II lines were 
measured by means of profile fits. The stellar parameters of the TOP and SGB group have not been corrected for helium 
diffusion, which would result in higher log g values (+0.05; see text). Below, the spectroscopic results are compared with 
the photometry (Stromgren and broad-band indices calibrated on the infrared-flux method 28 , see also Table 2) obtained 
with the Danish 1.54-m telescope on La Silla. In both cases, A log g is based on the magnitude difference in V. 



indicate for giant stars. However, the low signal- 
to-noise data did not allow firm conclusions. A 
different study 16 examined TOP stars and subgiant 
stars in NGC6397. It concentrated on elements 
showing suspected nuclear processing (oxygen, 
sodium, magnesium and aluminium). The analysis 
was based on standard models of stellar evolution 
and spectroscopic analysis (assuming equilibrium) 
and did not indicate significant differences in iron 
abundances. This result critically depends on the 
above-mentioned assumptions, which yield stellar 
parameters (effective temperatures) unsupported by 
photometry. The only previous study attempting a 
homogeneous analysis of red-giant and dwarf stars 17 
targeted M 13, a globular cluster a factor of three more 
metal-rich than NGC 6397. The highest-gravity object 
(a subgiant with log g = 3.8) indeed shows an iron 
abundance 41 % (0.23 in log abundance) below the 
average of the other 24 more evolved cluster members 
analysed. This was not interpreted as an indication 
of departures from a uniform cluster composition as 
caused by diffusion. 

To derive stellar parameters and elemental abundances 
from our observations, we employed well-established 
spectroscopic diagnostics (in one dimension) with a 
high level of modelling realism in the line formation 



(non-equilibrium where necessary). The profile of the 
Balmer line of hydrogen (Ha, at 656 nm) was used to 
estimate the effective temperatures (T e g) 18 ' 19 , while 
the ionization equilibrium of iron (treated in non- 
equilibrium) was used to determine the stellar surface 
gravities (logg) 20 . The mean stellar parameters for 
the four groups of stars are given in Table 1 . 
An independent check of the reliability of the 
stellar-parameter determination is obtained based on 
photometry, exploiting the luminosity difference of 
the stars. These data are also given in Table 1. As 
can be seen, the spectroscopic A T e s (TOP — RGB) 
(that is, T cff (TOP) - T eff (RGB)) is well-matched by 
both the Stromgren index (v — y) and the broad-band 
index (V — I), the photometry indicating a AT c ff only 
20-50 K (2-5%) lower. The spectroscopic Alogg 
(TOP -RGB) is only 0.05 below that determined 
from the photometry. This is excellent agreement 
for two fully independent methods and we thus 
consider the stellar-parameter differences to be well 
constrained. On the basis of these differences, relative 
abundance trends can be scrutinized. 
The best-determined abundance (in terms of number 
of spectral lines used) is that of iron. We find the 
abundance to differ by 45 % (0.16 in log abundance), 
that is, the TOP stars have the lowest abundance 
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Figure 1: Trends of iron and lithium as a function of the effective temperatures of the observed stars compared to the 
model predictions. The grey crosses are the individual measurements, while the bullets are the group averages. The solid 
lines are the predictions of the diffusion model, with the original abundance given by the dashed line. In b, the grey-shaded 
area around the dotted line indicates the la confidence interval of CMB + BBN 1 : log[e(Li)] = log (N-li/Nh) +12 = 2.64 ± 
0.03. In a, iron is treated in non-equilibrium 20 (non-LTE), while in b, the equilibrium (LTE) lithium abundances are plotted, 
because the combined effect of 3D and non-LTE corrections was found to be very small 29 . For iron, the error bars are the 
line-to-line scatter of Fe I and Fe II (propagated into the mean for the group averages), whereas for the absolute lithium 
abundances 0.10 is adopted. The 1<t confidence interval around the inferred primordial lithium abundance (log[e(Li)] = 
2.54 ± 0.10) is indicated by the light-grey area.We attribute the modelling shortcomings with respect to lithium in the 
bRGB and RGB stars to the known need for extra mixing 30 , which is not considered in the diffusion model. 



4 



Nature 442, 657-659 (10 August 2006) 



which successively rises towards the RGB stars 
(Fig. la). Propagating the line-to-line scatter for all 
iron lines of each star into the mean values given in 
Table 1 and further into the abundance difference be- 
tween TOP and RGB stars, the difference A log e (Fe) 
(TOP - RGB) is significant at the 3.2a level. Calcium 
and titanium also show trends (Fig. 3), but these are 
less pronounced. 

We have compared these abundance trends with 
various diffusion model predictions and found 
that a model with one particular value of the 
turbulent-mixing efficiency 10 is capable of fitting 
the observations of these heavy elements well. This 
model (model T6.0, with a parametric description 
of turbulent mixing using an isotropic turbulent 
diffusion coefficient 400 times larger than the atomic 
diffusion coefficient for helium at log T = 6.0 varying 
with density as p~ 3 (ref. 10), and computed for 
the metallicity of the RGB stars) clearly predicts a 
steeper trend for iron than for calcium or titanium, in 
agreement with the observations. The simultaneous 
element-specific reproduction of the observations thus 
lends strong support for the diffusion interpretation 
of these trends and constrains the level of turbulent 
mixing. 

We emphasize that it is not possible to remove 
all trends simultaneously by adjusting the stellar 
parameters: the neutral species Cal and Fel (both 
affected by non-equilibrium) require a T e g correction 
of roughly 100K and 200 K respectively, while the 
ionized species Ti II and Fe II (formed in equilibrium) 
would require a change in \ogg of 0.15 and 0.33, 
respectively. Considering the agreement between 
spectroscopy and photometry, these corrections to 
the stellar parameters are large and incompatible 
with one another. We also investigate the effect that 
the three-dimensional (3D) hydrodynamic nature of 
stellar atmospheres 21 might have on the observed 
trends. Using 3D TOP and RGB models available to 
us, we find that the trends inferred from weak Fe II 
and Till lines in one dimension require small 3D 
corrections only. For iron, the trend would even be 
slightly steeper. As the stellar-parameter differences 
are well determined, analyses based on 3D hydrody- 
namic models would reach much the same conclusion 
about the abundance trends. 

In the light of the identified diffusion signature, 
we should consider the structural effect that helium 
diffusion has on the atmosphere and, therefore, on the 
spectroscopic analysis. The T6.0 model predicts the 
He/H ratio in the TOP stars to be decreased by nearly 



50% from the original value, 40% in the SGB stars, 
whereas the original ratio is almost restored in the 
bRGB and RGB stars. Helium settling of this extent 
changes the mean molecular weight in the atmosphere 
appreciably, an effect which can be mapped as a shift 
in gravity 22 . The accompanying shift in log g amounts 
to +0.05 for both groups, bringing the spectroscopic 
A log g (TOP — RGB) into perfect agreement with 
the photometry-based value. This correction also 
improves the agreement of the stellar parameters with 
the isochrone (Fig. 2). 

Lithium is ionized easily and therefore behaves much 
like helium in terms of settling. As seen from Fig. lb, 
both the TOP stars and the SGB stars sample the 
lithium plateau. The observations support the T6.0 
model predictions, with a slight upturn towards cooler 
temperatures before the convection zone encompasses 
lithium-free layers and dilution sets in. A higher 
value of lithium among halo subgiants than among 
halo dwarf stars was recently reported 23 and qualita- 
tively interpreted as a signature of atomic diffusion 
counterbalanced by gravity-wave-induced mixing. 
In quantitative terms, the diffusion model we use 
predicts the original lithium abundance to be 78% 
(0.25 in log abundance) higher than the mean value 
measured in the TOP and SGB stars. Therefore, we 
infer a primordial lithium abundance which agrees 
with the cosmic microwave background plus Big 
Bang nucleosynthesis (CMB + BBN) value 1 within 
the mutual la error bars. Given the uncertainties 
associated with both values, the cosmological lithium 
discrepancy is thus largely removed. This finding 
restores confidence in standard BBN, quantitative 
spectroscopy and sophisticated stellar evolution 
models. 

In extrapolating to the primordial lithium abundance, 
we do not consider corrections related to the Galactic 
chemical evolution of lithium. Such corrections are 
model-dependent and rather uncertain. It is plausible 
that some fraction of the halo gas was processed 
through Population III stars 24 (lowering the lithium 
abundance) and that cosmic rays interacting with the 
interstellar medium have produced some lithium 6 . 
At present it is, however, unclear to which extent 
(and even in which direction) the stellar lithium 
abundances should be adjusted. 
While helium diffusion has been invoked to lower 
globular-cluster ages to values which meet the cosmo- 
logical age constraint, metal diffusion has rarely been 
considered. The models including atomic diffusion, 
radiative accelerations and turbulent mixing for all 
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elements will result in moderately different absolute 
and relative ages. Using similar models 25 , the analy- 
sis of M 92 has led to an age estimate compatible with 
Wilkinson Microwave Anisotropy Probe (WMAP 1 ) 
results. It is clear, however, that stellar-evolution 
models need to be developed further towards physical 
self-consistency with respect to all relevant processes, 
including hydrodynamics and rotation. In particular, 
the physical mechanism causing turbulent mixing 
should be identified. Promising work 26 points to- 
wards the importance of internal gravity waves in 
transporting angular momentum and mixing stellar 
interiors. 

Regarding unevolved metal-poor stars as tracers 
of Galactic cosmochemistry, heavy-element abun- 
dance ratios (for example, Mg/Fe) are less affected 
by diffusion than ratios with respect to hydrogen 
(for example, Fe/H). However, the effects are non- 
negligible and must be taken into account to reach 
the highest accuracy. Care should be taken especially 
when comparing abundances from giant and dwarf 
stars. In this respect, one exceptional pair is the two 
most metal-poor stars known 27 , whose abundance 
signatures will have to be reinvestigated in the light of 
diffusion. 
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Figure 2: Loci of the observed stars in the observational and physical parameter space. Top: Colour-magnitude 
diagram of NGC 6397 clearly showing the main sequence, the subgiant branch and the red-giant branch. The observed 
stars are indicated by the blue crosses. Bottom: Comparison of the spectroscopic stellar parameters (corrected for helium 
diffusion, see text) with a 13.5 Gyr isochrone constructed from the diffusion model T6.0. The agreement is satisfactory 
indicating that the absolute stellar parameters meet the cosmological age constraint. 
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Figure 3: Trends of calcium and titanium as a function of the effective temperatures of the observed stars compared 
to the model predictions. As weak lines were used, the analyses were done on the mean spectra for each group of stars 
using the mean stellar parameters given in Table 1. The calcium abundances (top) were determined in non-LTE from 
three neutral lines (at 612.2, 616.2 and 643. 9nm) which have line strengths between 30mA (Cal 612.2 in the TOP stars) 
and 92mA (Cal 616.2 in the RGB stars). The titanium abundances (bottom) were determined from the only available 
ionized line (Ti II 522.6 nm) which has line strengths between 22 mA (TOP) and 67 mA (RGB). For the relative abundance 
trends errors of 0.05 are adopted. The observed trends (markedly shallower than for iron) are in good agreement with 
the diffusion predictions (full-drawn line; the dashed line indicates the original abundances used in the model). Trends 
were observed for other elements (e.g. scandium, magnesium), but no model predictions exist for some (e.g. scandium) 
and others (e.g. magnesium, aluminium) suffer from suspected nuclear processing in globular clusters 16 and are thus not 
suitable for isolating the diffusion signature. 
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Spectroscopic versus photometric effective temperatures of the four groups of stars observed 



Group 


No. of stars 


T cff (spec) (K) 


T cS (v - y) (K) 


T cS (V - /) (K) 


TOP 


5 


6254 ± 9 


6240 ± 21 


6133 ± 31 


SGB 


2 


5805 ± 7 


5824 ± 6 


5688 ± 8 


bRGB 


5 


5456 ± 57 


5408 ± 37 


5318 ±29 


RGB 


6 


5130 ± 19 


5132 ±22 


5063 ± 20 



Table 2: The photometric indices are calibrated on the infrared-flux method . The agreement between spectroscopy 
and (v — y) is excellent, whereas (V — I) is offset towards cooler effective temperatures by roughly 100 K. However, 
for the identification of the diffusion signature the effective-temperature difference is the most relevant quantity. The 
effective-temperature difference between TOP and RGB stars in (V — I) is lower than in the other two cases (by roughly 
50 K). Because of this, a slightly larger abundance difference would result when referring to the (V — I) temperature 
scale. The conclusion is that three prominent effective-temperature scales (Ha spectroscopy, Stromgren and broad-band 
photometry) reveal significant abundance differences between the TOP and RGB stars of NGC 6397. 



